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. Numerical Simulations of the random phase sine-Gordon model suffer from strong finite size 

effects preventing the non-Gaussian log'^ r component of the spatial correlator from following the 



^ \ universal infinite volume prediction. We show that a finite size prediction based on perturbative 

, Renormalisation Group (RG) arguments agrees well with new high precision simulations for small 

' coupling and close to the critical temperature. 
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I. INTRODUCTION 

The random phase sine-Gordon model plays a role in understanding a variety of phenomena. The model has been 
hysically interpreted in terms of roughening of crystalline surfaces [1] , pinning of flux lines in superconducting film 



•ated a variety oi 

i-rj I predictions ll|, |5|. | 10| . not all in agreement, and continues to be a scene of activity [H [11 [11 [H. Numerical 
: simulations [131114 lla. [H [13, [H, 111, IH [111 might be expected to play a decisive role in this situation, but this 



I [1, [3) [E. [3 charge density waves ^, JJ, and other systems. The model has generated a variety of theoretical 

O has not in fact been the case. Over the past few years, some of the differences between simulations and theory have 

^ O ^ , been resolved as a result of better simulations and analysis [IB, [l3| and corrections to the precise values of theoretical 
predictions [H, [IB], but there still remain discrepancies. 

] The primary discrepancy between theory and simulation concerns the presence of a log^ r component of the spatial 

^ ' propagator computed at distance r in the high temperature phase (since the high temperature phase is Gaussian, the 

On [ spatial correlator should be simply logr). In the renormalisation group (RG) approach, this component is expected 

£0 • below the critical temperature at length scales beyond a characteristic length that diverges as the critical temperature 

^ ^ '■ is approached from below. However, in finite lattice simulations, this component is observed (isj above the expected 

I . critical temperature, depending on the strength of the disorder, and leads to uncertainty in actually identifying the 

ly-j ' location of the critical point . Although t he p ossibility that the critical temperature could depend on the strength 

' of the disorder was raised in theoretical work [2J|j the accepted theoretical view is that this does not occur. The reason 
for this is a symmetry of the Hamiltonian, noted by Cardy and Ostlund (25| and recognised as a particular example 

■ of the so-called statistical tilt symmetry [2a] . Another discrepancy concerns the strength of the log^ r component in 
^ i: the low temperature phase. The RG prediction for this strength turns out to be universal, yet it does not match the 

I ■ value observed in simulations on achievable lattice sizes. Figure ^ illustrates the discrepancies between the infinite 
volume RG prediction and data from simulations of several years ago [isj . These phenomena have been observed in 
other numerical studies [T3| . 

^ , In this paper we use perturbative RG arguments to estimate the size of the log^ r component for small couplings in 
IL" ■ the vicinity of the critical temperature bearing in mind the finite size of the lattices used in simulation. We confront 
. 5^ [ these estimates with data from new high precision simulations made within the restricted parameter range where 

■ perturbation theory can be relied upon. Our estimates do display a log^ r behaviour above the critical temperature 
^ ' and match the strength below the critical temperature therefore resolving the discrepancy. 

To use the RG theory to make predictions on the lattice, we need to compute the coefficients appearing in the RG 
equations using a lattice regulator. The usual treatment is via a Coulomb gas |25l] which involves a completely different 
type of regulator. We have obtained RG coefficients for the lattice regulator by adapting Kogut's original approach 
of momentum shell integration [23|. These coefficients are substantially smaller than the conventional Coulomb gas 
ones, and allow quantitative agreement with the simulation data. Nonetheless, the RG theory is only valid for small 
couplings, and only in this region do we obtain good agreement. 
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FIG. 1: Numerical results for the coefficient 62 (see the text) of the log^ r term in the correlator on a Q4^ lattice using data 
originally discussed in [itI. [T8||. Bottom to top: A = 0.5,2.0 and the discrete height model (formally equivalent to the sine 
Gordon model with A = 00). The vertical line indicates the position of the infinite volume critical temperature at Tc = 2/n, 
and the horizontal line shows 62 ~ 0. The wider line marks the (universal) infinite volume RG prediction (&2 ~ 2r^). 



We commence with a description of the model and simulation data we wish to explain but relegate technical details 
of numerical work to appendix A. Then in section III, we briefly review the RG for this model and refer to appendix 
B for the full computation of RG coefficients appropriate to the lattice regulator. Section IV discusses our technique 
for estimating the finite size correlator and presents results in the form of a figure confronting the prediction against 
the new simulation data. The conclusion critically discusses our treatment and since this work does not constitute a 
full finite size theory, it also discusses some of the issues involved in such a theory. 



II. SIMULATIONS 



The random phase sine-Gordon model has Hamiltonian, 

^ = II (<^» - )' - ^ II cos 27:{(b, ~ m) . (1) 

Depending on the physical system, the continuous field </) may be interpreted in different ways. In the case of a surface, 
(/> is the height of the interface so the field is not periodic and vortices are excluded. The disorder is contained in the 
term r\ which is characterised by some probability distribution. Usually a flat distribution for the quenched disorder 
in < 77 < 1, is assumed, but other distributions have also been simulated [l3| . 

In order to clarify the problems that arise in simulations, let us start by recalling the results of our earlier numerical 
study [l^ of this model which form the basis of figure Signals of a transition were clear in both the static 

and dynamic properties [38| . Although the static correlation function had the expected logr behaviour for high 
temperature, it was clear that another component was present at low temperatures. The results were not sufficiently 
sensitive to determine the functional form of this additional component, but a log^r term as predicted by the RG 
studies allowed a good fit to the data. To be precise, the simulated spatial correlator was fitted with the function 
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(using ((• • •)) to denote average over the thermal noise and (• • •) as the average over the quenched disorder): 



G{r) = {{4>r - M ) = b,PL{r) + b^Pl{r) (2) 
hi and &2 are the fit parameters and PL{f) is the Gaussian correlator on a lattice of size L, 

L-l L-l 1 _ pQg('2£rrii ^ 1 

m— ln2— ^ L ' ^ L ^ 



and the symbol ~ holds for large lattices, L ^ 1 and length scales, r, that are large but not approaching L. 7^ is 
the Euler-Mascheroni constant. This is actually the correlator along a lattice axis and is the appropriate quantity to 
compare with the correlators obtained from simulation. 

Instead of relying on our old data to make comparison with our predictions we have performed new high precision 
simulations. These new simulations at small coupling, A = 0.5, remain on the same 64^ lattice but have better 
thermalization, an improved random number generator and higher statistics. Details of the simulations and comparison 
with the old data are discussed in Appendix A along with additional data exploring other values of A. Data for the 
coefficient 62 of the log^ r term based on new simulations are presented in figure ^ later in the paper. 



III. RENORMALISATION GROUP 



A convenient basis for the RG studies is the continuum replicated version of the Hamiltonian ([T]) following Cardy 
and Ostlund [H: 

13H= fd^x^Yl Kc.[)d(t>o.d<j>(i - ^ ^ cos 2^(0„ - ^p) . (4) 

Greek letters denote the n replica indices and we have made g dimensionless by explicitly writing the UV cutoff a. 

The kinetic term of the replicated Hamiltonian is parameterised with separate coefficients for the on and off-diagonal 
pieces, Kaf^ = KSaf^ + {K — K){\ — Sap), since the off-diagonal piece appears after renormalisation even if not present 
initially. The bare Hamiltonian only has a diagonal term, so the bare parameters are: K{0) = K{0) = 2/3. The bare 
value of g is related to the parameters of the simulated Hamiltonian by: g{0) = (/3A/2)^. 

RG equations may be obtained from ^ [H, [2^ [3l|, [3^, HI], though the same results are obtained without 
replicas [5|. In either case an additional term is generated: in the replica calculation it appears as the off-diagonal 
part of Ka/s, whereas in the other approach it appears as the coefficient of a new term in the Hamiltonian, linear 
in the derivative of <j). These same equations have been derived in a variety of ways; a Coulomb gas approach [25[, 
conventional field theory [sij , conformal field theory (33| , with the help of the exact renormalisation group [34j| and in 
appendix A, we ourselves use momentum shell integration. Irrespective of the method of derivation, the parameters 
change with the RG scale, I, in the following way: 

dg^ 

dl 
dK 

'df 
dk 
'df 

Where the reduced temperature r = (1 — Kc/K) 
Kf. — tt). These perturbative RG equations have been obtained in a double expansion: small coupling g and small 
reduced temperature r (3lj | . They are valid for small couplings in the region near the critical point where higher order 
terms can be ignored. 

The asymptotic form as Z — s- 00 of the solution of these equations depends on the phase. In the high temperature 
region where r is negative, the fixed point value of g{l) is zero and K{1) goes to a constant. In the low temperature 
region g{l) becomes g* = "^t/C, and the coupling K{1) behaves as K{1) — — 4(A/C^)t^L The consequences of these 
solutions are fully explored in the literature and we do not repeat the analysis here. 

The RG coefficients A, C are dimensionless, non-universal constants of order unity. They depend on the temperature 
but various workers [1, [3l| noted that the amphtude ratio A/{KC)^ appearing in the asymptotic form oi K{1) below 



= 2Tg~Cg\ (5) 
= ~Ag\ (6) 
= 0. (7) 
(1 — T/Tc). The critical temperature is at Tc = 2/7r (hence, 
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Tc, is universal at the critical point. The correct value for this ratio was clarified by Carpentier and Le Doussal |23| 
(see also reference [l^)- In our conventions (note that we keep a factor of 2tt inside the cosine term), this ratio is: 

A R T 1 T 



{KC f (27r)2 Tc (4^) ' 

R is the universal ratio defined in [23| and carefully shown to take the value tt. 

The non-universal values of A and C themselves are necessary to compare with simulation. The regulator natural 
in the Coulomb gas approach used by Cardy and Ostlund [IBl gives values C = 47r and A = A-k'^. These values are 
rather high and we should properly use values from a lattice regulator. A full perturbative RG computation on the 
lattice becomes difficult so we have resorted to the following argument. 

We have used momentum shell integration following Kogut '^t'I to rederive the RG equations. This technique 
needs a modification first pointed out in 35] to take account of the correct operator product expansion, however all 
computations can be consistently expressed in terms of Gaussian integrals without any need to be concerned with 
the subtleties of two dimensional massless canonical field theory. Details are given in appendix B, and this technique 
yields an expression in terms of the Gaussian propagator < (j)a{r)<j)p{^) >= K~pGQ{r): 



C = / A^^e^l«"(«)-'^''(")1 = 47re*(°°) . (9) 

K J dK 

The final expression on the right is evaluated at the critical point, is identical that obtained by a different procedure 
in [23I I and gives the correct value of the universal ratio. $ depends on the regulation scheme and is defined in terms 



of the subtracted propagator: 



[Go(0 - Go(0)] - log(rA) + ^<^{rK) . (10) 



The momentum shell integration technique requires a regulator that is sufficiently smooth to make the term '^'^^j^^ 
short range and thus render the expression for C well defined. The lattice propagator has this property, and we can 
identify $(00) using the asymptotic form of the lattice regulator ([3]) [3^ : 

[Go(r) - Go(0)]i,,,i,, ^ -;llog(rA2V2e^-) . (11) 



We therefore use the value of C given by 



C = -e-^^^ . (12) 



Recalling the value of the Euler-Mascheroni constant 7^ ~ 0.577..., we find this value of C to be substantially smaller 
than the Coulomb gas value. 



IV. FINITE SIZE 



In order to compare the RG theory with simulation, we need the correlation function on a finite geometry. The 
usual derivation of the RG correlator in infinite volume was given by Toner and Vicenzo [l| and we shall proceed in 
a similar way without taking asymptotic values for the RG scale. Numerically it is better to work in the real space, 
but the fields (f) do not transform simply under the RG, so we consider the composite fields e*^"^ that have better 
properties [Slj. General RG arguments allow us to write an equation for the correlation function Cg(r) of these vertex 
operators (see for example reference [3l|): 

C,(r,i^,if(O),g(O)) = (exp(jg(0(r)-0(O)))=exp(- / 7,(g(x))— ) C,(l, ii^, i^(logr), g(log r)) . (13) 



l/r 



X 



Where ^q{g) is twice the (RG) dimension of the vertex operator e-xjpiiqcfiir) . The leading terms in the perturbative 
expansion of this quantity are: 
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Determining the value of 72, twice the anomalous dimension, in our RG scheme is not straightforward (it requires 
the OPE of three vertex operators) and neither has it been computed in other renormalisation schemes. We have 
therefore only kept the leading term in 7^ when expanding equation (|13p in q. We expand to order and write 

Cq{r) = 1 — ^G(r) + 0{q'*'), to obtain an RG equation for G{r): 

f^°s{r) 1 / K(l)\ 
G{r,K,Ki0),gm = G{l,K,Ki\ogr),g^0) + 2j^ (^2 - -4^ j . (15) 

This expression amounts to using the zeroth order formulae improved with the running coupling constants computed to 
second order in perturbation theory and is the same procedure usually followed. We will show that this approximation 
provides us with the same qualitative picture obtained in numerical simulations and a good agreement with the 
numerical data obtained at small coupling constant. Besides omitting the order term from the renormalisation 
factor 7q, we have consistently terminated the perturbative expansion of G at zeroth order. In fact we have obtained 
an expression for the order g^ contribution to G, but the resulting formula is not amenable to numerical calculation 
on suitable size lattices [S^. We remark that even if we were able to include these order g^ terms, it would only be a 
step towards greater precision since the complete model on the lattice also requires irrelevant operators (see appendix 
C). 

On the lattice we replace logr by its lattice version equation ([3]), log(r) 2'K{PL{r) — Pl(1)), so our final equation 
for the correlator is: 

As is the case for the momentum space version, this equation depends on the running coupling K{1) and consequently 
the behaviour of the correlator depends on the phase. The standard, infinite volume, analysis is based on using the 
asymptotic solution for K{1). In the high temperature phase where K{oo) is constant this yields a logr correlator 
with the correct coefficient {T/2t:) even at temperatures beyond the region of validity of the perturbative result. At 
low temperature the asymptotic solution K{1) — 4(A/C'^)r^Z leads to the log^ r behaviour [l|. 

G{r,K,K{0),g{0)) = 2r2p|(r) + ... (17) 

This leading term has universal coefficient, which takes a value 62 = 2t^ (see equation ([2|)). The sub- leading logr 
term, (the dots in the equation above) can be obtained from corrections to the asymptotic value of K{1), though its 
precise value is dependent on the regulation scheme. Comparison of the two terms yields a length scale that separates 
the logr from the log^ r behaviour. Since the coefhcient of the log does not have any r factor, this crossover scale 
diverges as e'^°'^''*/'^ as Tc is approached from below. This is the problem identified in the Introduction. 
To make predictions for finite size, we use the full solution of the RG equation for K{1) psj : 

K{1) = K{0) Dr (log (1 + x{e'^' - 1)) - x(l - x) J^^^j^iTT^) ■ (18) 

Where the parameter x = 9(0)G/2t has the same sign as r. D — 2A/C'^ = 1/T, is universal, as discussed in the 
previous section. This is an exact solution of the perturbative RG equations valid in the region of small g{l), namely 
for small A near the transition. By using this expression and numerically solving the real space RG equation (i.e. 
computing numerically the integral which appears in the right of equation (|16p ). we are able to compute the finite size 
lattice correlator G(r). The resulting correlator could be compared directly against simulation data for this quantity, 
but it is clearer to present the comparison in the same way as we did in our earlier work by fitting the predicted 
correlator with a form ([2]). 

To put this program in practice we use the value of G = ne"'^'"^ /2 that we argued was appropriate for lattice 
regulation, and we make the fit to the finite size correlator over the range r = 1, 32 since both our original and new 
simulations are on a 64^ lattice. The resulting predictions are shown in figure These predictions have the same 
qualitative form as the simulation data shown in figure ([1]). Namely, we observe the same apparent shift in critical 
temperature with A or appearance of log^ above Tc and the same dependence of this shift, which increases with 
increasing A. 

For a more precise test of agreement, figure ([3]) expands the previous figure near the critical temperature to 
confront the new high precision A = 0.5 data against the predictions. The features identified in the Introduction 
are clearly apparent in the data in figure ([3]) : the continued presence of the log^ r coefficient above the expected 
Tc = 2/7r = 0.6366, and the disagreement with the universal infinite volume prediction for this coefficient below 
Tc- The similarity between the finite size prediction and data is striking. Within the temperature range shown, the 
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FIG. 2: Theoretical estimate for the coefficient 62 of the log^ r term in the correlator as a function of the temperature as 
determined by fitting RG predictions for a lattice size 64^ using the method described in the text. We show continuous lines for 
five A values (bottom to top in the plot): 0.1, 0.2, 0.5, 0.7 and 1.0. The vertical line marks the position of the infinite volume 
critical temperature. The wider line marks the (universal) infinite volume RG prediction (62 = 2r^). 



predictions lie within twice the error bound of the simulations. We can quantify the extent of agreement by computing 
between the numerical data and our theoretical prediction. We obtain = 5-13 using the temperatures in the 
interval [0.5, 0.7] , the number of degrees of freedom (d.o.f.) being 5 (notice that we fit the numerical data to a function 
that we have obtained without resort to these data). The probability to have a larger than = 5.13 for 5 d.o.f. 
is 40%, which is larger than the 5% limit typically used in the literature to discriminate the validity of hypotheses 
[S^]- If the interval is enlarged further below the critical temperature to [0.45,0.7] we obtain — 7.85 with 6 d.o.f. 
corresponding to a probability of 25% (40| . This agreement can also be seen by directly comparing the prediction for 
the correlator against the result from simulations. At T ~ 0.55 the agreement is very good, with all predicted points 
lying well within the small errors of the data. 

In Appendix A, by way of motivating our choice of parameters for the new simulations, we discuss some exploratory 
simulations at higher values of A. These suggest that for larger couplings as we move away from the perturbative 
regime, the agreement will not be so good. Indeed, in figure ([T]) the A = 2.0 and A = 00 data are not far apart, 
suggesting that A = 2.0 is well into the non-perturbative regime. For these reasons we did not perform any further 
new simulations. Nonetheless, the predictions in figure ^ retain the qualitative features observed in the original 
simulations, figure ((T|). 

Although our simulations were only at a particular lattice size, it is illuminating to consider how the results of 
the procedure above depend on the lattice size. We find that above Tc, as the size is increased, the coefficient 62 
decreases, but the dependence is extremely weak as was found in numerical simulations (in references [T7l.[l^. L = 64 
and L = 128 lattices yielded 62 curves with overlapping error bars in the low temperature phase). By following our 
finite size procedure for a range of lattice sizes up to 512, we find a good fit to a scaling form. At T = 0.7 the 
coefficient 62 goes to zero according to l/L"-^^ for both A = 0.5 and 1.0. This exponent takes a value close to 4t since 
at this temperature r = —0.101. This result corresponds to a finite size dependence b2{L) oc g{L)'^, that might be 
expected since g enters the RG equation for K as g^. 
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FIG. 3: Comparison of theoretical estimates for the coefficient 62 of the log'^ r term in the correlator against new high precision 
simulations at A = 0.5 in the vicinity of the critical temperature. Axes and marked lines have the same significance as in 
earlier figures. The continuous line is the theoretical estimate for A — 0.5 and the points are from the new A = 0.5 numerical 
simulation (L = 64). 



V. CONCLUSION 



Our main result lies in the similarity shown in figure ([3]) between new simulations and the predictions based on the 
RG estimate according to our finite size treatment. Until now, the numerical simulation results have been a puzzle 
since it appears that either there is a log^ r component to the correlator above the critical temperature, or that the 
location of the critical temperature depends on the disorder. Neither of these interpretations are consistent with 
the conventional RG treatment. In this paper we have shown that one can qualitatively understand the numerical 
observations in the framework of the RG without the need for any additional ingredient just by taking account of 
strong finite size effects. Above and near the critical temperature, the disorder strength renormalises to zero and the 
fixed point is Gaussian. However the dependence on lattice size is extremely slow, and this behaviour induces the 
log^ r term above the critical temperature. 

The statistically significant level of agreement relies on several factors. Firstly, the use of new high precision 
numerical simulations that are better thermalized than the old data. Secondly, the agreement only occurs in the regime 
where perturbative RG is valid, namely for small disorder strength (A < 0.5) and near the transition 0.5 < T < 0.7. 
Finally the level of quantitative agreement is crucially dependent on the use of RG coefficients appropriate to the 
lattice regulator rather than the traditional Coulomb gas values. 

Our treatment does not constitute a full finite size theory, and there are several factors we have neglected. The first 
of these is the perturbative correction from g{l) to the equation (jlSp . Although the leading perturbative contributions 
to the correlator can be written down, their numerical evaluation would require additional approximations that 
would undermine attempts to make detailed quantitative comparisons. Moreover, the renormalisation factor from 
the anomalous dimension of the vertex operator has not been computed. A proper treatment of these perturbative 
corrections would require a full lattice based perturbative computation. Another factor is the contribution from 
irrelevant operators. Irrelevant operators can arise from the higher cosine modes after integrating out the disorder 
and also from the discretisation. A brief discussion of these operators is given in appendix C. A final factor, the use 
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of an RG coefficient appropriate to the lattice regulator has been included in our treatment, and indeed without this 
we would not achieve anything like the quantitative agreement shown in figure ([5]). 

In our earlier numerical work [l^ we also investigated the dynamics. We found that the dynamical critical exponent 
z departs from its Gaussian value above the critical temperature, and that the value of this term depends on the 
disorder strength (see also 0]). We expect that this effect can also be explained in terms of the same strong finite 
size effects. 
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APPENDIX A: SIMULATION DETAILS 

Here we discuss technical details of both the old and new simulations and present additional numerical results. All 
the data under discussion has been placed in the Journal repository. 

1. Original Simulations (dataset A) 

Our original numerical simulations [isll were performed in 1995 on an APE supercomputer. Working on a 64^ 
lattice, two values of the coupling were considered, A = 0.5 and A — 2.0 and 128 samples were used for each data 
point. Figure ([T]) also shows results for the discrete model [13, [HI, but in view of the large coupling this model 
corresponds to, we do not consider it in the present work. 

The data were based on a single annealing run (henceforth called dataset A), with 10^ thermalization and 2.10^^ 
measurement sweeps at each temperature (separated by AT — 0.05). The random number generator used in this 
work [28| has since been criticised [1^ . 

2. New Simulations (datasets B,C) 

To reduce the size of the error bars on dataset A so as to make quantitative comparisons we performed a number of 
high precision simulations around the transition temperature. Since the APE machine used for the original simulations 
is no longer available, we wrote new code running on XEON PC's. These new simulations remain on the same 
64^ lattice but have better thermalizatioUj higher statistics and an improved random number generator (a 64-bit 
congruential generator based on reference ^30]). Both dataset B and C are at A = 0.5. 

Dataset B checks the new code by using the same parameters and the original annealing schedule of dataset A and 
closely reproduces those results. Error bars generally overlap, though they are slightly apart at T=0.65,0.7. A ^ 
evaluation gives a confidence limit of 15% that the data arise from the same distribution. The new data (B) tends to 
give slightly higher signal than the old data (A) and we would ascribe this difference, if significant, to the improved 
random number generator. We therefore have confidence in the new code. 

For improved thermalization, dataset C, the annealing schedule has been lengthened and consists of separate runs 
of 500 samples with 10^ thermalization sweeps at temperature intervals of AT = 0.1. The final measurements are then 
taken over 10^ sweeps. Each data-point (C) lies within two sigma of the other data, however, there is a systematic 
effect that always increases the new value of hi relative to either A or B. Thermalization is the main factor in the 
difference, but we cannot discount an effect also from the random number generator. 

We have also performed runs to observe the approach to thermalization at T = 0.45 by tracking the value of the 
correlator at maximum distance (32, half the lattice size), as the measurements progress. This study indicates that a 
thermalization of 10® sweeps is necessary and sufficient (to an accuracy corresponding to the number of samples we 
use). We are therefore confident in the new data C and this is what we use in the main text to make quantitative 
comparisons with the RG predictions. 
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A 






<x 


0.085(5) 


- 


2.0 


0.067(6) 


0.268 


0.7 


0.035(2) 


0.048 


0.6 


0.031(2) 


0.033 


0.5 


0.018(2) 


0.020 


0.4 


0.012(1) 


0.0105 


0.3 


0.006(1) 


0.0041 



TABLE L Data for 62 just below Tc at T = 0.50. Simulation data (original for A = 00, 2.0 other values from new simulations) 
denoted by fe""™ compared with predictions denoted by 6*. 



A 


I num 


7 th 
02 


00 


0.025(3) 




2.0 


0.018(4) 


0.108 


1.0 


0.0050(15) 


0.020 


0.7 


0.0046(17) 


0.007 


0.5 


0.0034(9) 


0.002 



TABLE II: Data for 62 just above Tc at T = 0.65. Simulation data (original for A = 00, 2.0 other values from new simulations) 
denoted by fef"™, compared with predictions denoted by b^. 

3. New Exploratory Simulations 

Before deciding the coupling strength to use for run C, we made a short set of new simulations to explore new 
values of A at two fixed temperatures: T = 0.55 just below Tc and T = 0.65 just above Tc- 

Table U shows results for a variety of A just below the critical temperature at T = 0.50 (128 samples in each case). 
At this temperature we are able to make accurate measurements at small A and have included the data from both 
new and old simulations. We find that for A < 0.6 the predictions agree well with the results of the simulations, but 
at A = 0.7 the value is not statistically compatible. 

Table HIl shows similar data, based on 500 samples, for a variety of A just above Tc at T = 0.65. For A > 0.5 the 
system shows a clear signal of log^ r, but we have been unable to obtain statistically useful results for smaller A due 
to the number of samples required. We find agreement with predictions at A = 0.5, though not at A = 0.7, matching 
the range of agreement just above Tc- 

On the basis of the fixed temperature exploratory simulations it appears that A = 0.5 is within the range of 
small coupling where the RG computations are accurate. We therefore collected dataset C at fixed A = 0.5 over the 
temperature range [0.4,0.7]. 

APPENDIX B: RG COEFFICIENTS BY MOMENTUM SHELL INTEGRATION 

In this appendix we use momentum shell integration technique to derive expressions for the RG coefficients that 
can be used with a lattice regulator. Although we end up with formulae equivalent to those obtained using other 
methods in [23| , it is only by going through the full derivation that we are able to check that the expressions remain 
valid with a lattice regulator. The momentum shell technique is direct and intuitive and is an interesting calculation 
in its own right. The technique is explained by Kogut for the quite similar (ordered) sine-Gordon model in [27*1. Since 
we follow this reference quite closely, in this appendix we only give an outline of the computation with particular 
details of differences from the ordered sine-Gordon case. 



1. Notation 

It is convenient to introduce some notation for the local operator: 

Oc,[3{x) = e2^H0=W-0/3W) . (Bl) 
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Also for quantities that appear after Gaussian integration. 

Where G{x) is the full propagator of the fields and Gh{x) is the propagator for the h fields that have restricted 
support only in the momentum shell (A — dA < p < K). 

G,{x)^dA—. (B4) 

The replica indices only appear through the inverse kinetic matrix and only take distinct values on or off the diagonal. 
It is therefore useful to write: 

A{x) = Aaffix) a = fi (diagonal) , (B5) 
A(x) = Aafj{x) a 7^ /? (off - diagonal) . (B6) 

And similar expressions for B and B. 

Then, using the explicit expression for A'^^, we find a simple form for the combination: 

A{x)A-^{x) ^ e.-'^^"-GK(x)lk _ (37) 

In which K does not appear. There is an analogous expression for the i3's. 

We also define coefficients 6a/3-y5(^) to take account of the connected expectations at second order, given by: 

6a/375(.T) = Al^{x)Al,{x)All{x)A^,^^{x) - 1 . (B8) 
2. Expectations 

The change in action as fields (denoted as K) in the momentum shell are integrated out is evaluated perturbatively 
in g. At first order the Gaussian integration yields the intuitive contractions: 

(cos2^(^„(a;) - <^^(.t) + h^(x) - /i^x))) - \ (o^p{x){^^'^''-^^'^~''^^-'^^) + c.c.) 

= yl2(0)i2(o)cos27r(</.c«(a;)- 0/3(2;)). (B9) 

The second order term is not much harder. In this and all subsequent computations we do not worry about operator 
ordering or normal ordering and always follow the consistent prescription of the path integral. 

(cos27r ((?!)a(a;) - 4>!i{x) + ha{x) - hp{x)) cos27r (0^(y) - ^^(j/) + h^{y) - /i5(y))) connected 

= ^A\0)A-HO) ( W(0Oa/3(a:)O^5(y) + c.c.) . (BIO) 

The cosine operators are at different spatial points x and y and we have written the difference: £^ = x — y. 

The coefficients bap-ysiO take account of the connected form of the expectation and were defined above. For a 
smooth cutoff, we expect the b{£^ys to be short range and allow us to use an operator product expansion (OPE) for 
Oaf3{x)0~fsiy)- Indeed, the momentum shell approach requires a regulator function that allows this step to proceed. 
The OPE was not needed in Kogut's work, but as pointed out in [35l|, is necessary in general. 

3. Operator Product Expansion 

The relevant terms in the OPE are: 

Oap{x)0^s{y) ^ ai{^)5asSi3y 

We compute the coefficients ai and 02 consistently by using exactly the same Gaussian integration techniques we 
have employed throughout the computation. 

ai(e) = -Ti^B\Q)B-\Q)eB-\OB\0. (B12) 
a2(0 = B\Q)B-\Q)B-^{S,)B\0- (B13) 
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4. RG coefficients 



After rescaling the action, we identify the changes to couphngs necessary to preserve the physics, and obtain 
expressions for the RG coefficients. 



C 



(n-2) / d^^a^iO A-\OA\0-^ 



K 

k 



dA 



-e K 



-[G(i)-Gm 



d^^A 



dA 



■[G«)-G(0)] 



(B14) 
(B15) 



Notice that these expressions are in terms of a regulated propagator and that although the raw propagator G{x) is 

Fra-red divergeni 

Having checked 
a lattice regulator 



infra-red divergent, the formulae appear in terms of well defined propagators: ^^^^ and G(^) — G(0). 

Having checked that ^^^^jp- is short range for the lattice propagator, we can be confident in using the approach with 



5. Smootii cutoff and Universality 



Recognising that G{r) — G(0) is a function of the combination rA and that dG{0)/dA = l/27rA; we can write the 
dG{x)/dA part of the RG coefficients as: 



dG{x) 
dA 



d 

dA 



[G{x) - G(0)] 



1 



27rA 



(B16) 



Then re-expressing the d/dA as a a;.V, the derivative can then be taken to act on the exponential to obtain (in the 
spherically symmetric case): 



K 

A = l^A^ 
K 



d^e 



. K ^d 



e K 



■[G(«)-G(0)] 



. K , d 
'^-Jdl 



g^[G(«)-G(0)] 



(B17) 
(B18) 



Now, without taking limits of any kind, write the propagator in the exponent in the following form for any value of r. 

(B19) 



G(r) - G(0) = log(rA) + -^<f (rA) , 
27r 47r 



This is a definition of $ which contains all the regulator dependence. For small r, this will appear to be an unhelpful 
decomposition, since $(r) must contain terms to cancel the singularity in the logarithm. For large r, the logarithm is 
the leading term and all remaining regulator dependence comes from $(oo). This $(00) is the next to leading term 
of the asymptotic expansion of the propagator that is merely a constant and is often incorporated into the logarithm. 
Evaluating A and C at the critical point {Kc = tt): 



G = 47r 
A = 47r3 



dx- 



dx 



■I'(oo) 



dx 



(B20) 
(B21) 



Our expressions for the RG coefficients confirm the universal nature of the ratio, since the regulator dependent 
pieces contained in $(00) cancel in the ratio. Moreover the ratio takes the correct [23| value, R = tt. 



APPENDIX C: IRRELEVANT OPERATORS 



In the theoretical study of the disordered sine Gordon model one usually picks up only the relevant operators (in 
the RG sense) . However the irrelevant ones can induce scaling corrections which modify the behaviour of the system. 
The goal of this section is to identify such operators and to show up their influence in the scaling of the numerical 
observables. 
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We can identify two main sources of irrelevant operators. The first one is the discrete nature of the model which 
leads to: 

0(r + ai") = cPir) + a5^<A + ^a^d^^cp + ^a^d^^ + 0{a^) , (CI) 

where a'^ is the lattice vector with modulus a. Hence, in addition to the original kinetic term in the Hamiltonian 
(9^^)^), the following terms appear, 



92 



J d'^x {dl^f+g^ j (fx {dlct>f+g^ J d^x (a^0)(9^^)+ff5 J d^x {d^ml<l^)+9^ J d^x {dlcj,){dlcj,)+0{a'') . (C2) 



One can compute the RG equations to the leading order for each coupling by simply computing their (momentum) 
dimensions using the Gaussian Hamiltonian: 

^ = (dim 5,) 3, + ..., (03) 

with dim5i2 = —d, dim52 = —d, dimgs = —d — 2, dim 54 = — rf+ 1, 6xmg^ = —d, dim^g = —d— 1. Their value at the 
lattice scale L is given by: 

gi{L) ~ . (04) 

The second source of irrelevant operators is the disorder average. In this paper we consider a flat distribution for 

the disorder. In addition to the two replica term, the full disorder average of the model induces couplings with an 
arbitrary number of replicas. The next to leading terms, with four replica interactions, are: 

Wl ^ y d'^X C0s2-K{(j)a + (f>b - (f>c - <t>d) + U2'^ J d'^X C0s2'K{(f)a - (f>b) cos2'k{(1)c - (fyd) ■ (C5) 

ahcd abed 

By computing their dimensions with the Gaussian Hamiltonian, one can write the leading order RG equations for u\ 
and U2'- 

du ' 

-^ = 2{2T-l)ui + .... (06) 

Therefore both of them are irrelevant (from) well below the critical temperature. At the critical temperature one 
should expect scaling corrections proportional to 

Ui{L) ~ . 

A similar computation shows that the six replica terms would induce corrections proportional to L^^. In general, 
an n-replica interaction induces corrections proportional to Olearly, all these scaling corrections, are much 

weaker than the leading (two replica) one, g{L) o± i"^"^. 
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